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Abstract 

04 ' The system of extended objects interplaying with a black hole describes or mimics various grav- 

■ itational phenomena. In this brief paper, we report the results of stability analysis of codimension- 

one Dirac-Nambu-Goto branes rest at the equatorial plane of d-dimensional spherical black holes, 
including the Schwarzschild and Schwarzschild-(anti-)de Sitter black holes. For the Schwarzschild 
and Schwarzschild-anti-de Sitter backgrounds the stability of branes is shown analytically by means 
of a deformation technique. In contrast, for the Schwarzschild-de Sitter background we demonstrate 
, with the help of numerics that the brane is unstable (only) against the s-wave sector of perturba- 



OO 

o 



o 



X 



tions. 



1 Introduction 

The system of an extended object such as a string and a brane interplaying with a black hole appears 
in various contexts of gravitational physics. As a first stage of this direction, one may customarily 
neglect the thickness and the backreaction of brane into the geometry, in which case the motion of 
brane is governed by Dirac-Nambu-Goto action [T]. In spite of this test brane approximation, the 
system, which will be referred to as a brane-black-hole system, is able to provide a rich variety of 
physically intriguing results. Prime examples are the cosmological domain wall interacting with a 
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primordial black hole, a domain wall pierced by a black hole [2], and a probe Dg-brane interacting 
with background Dp-branes, and so on. The final topic quoted above has attracted much interest 
from the viewpoint of the gauge/gravity duality, according to which Dp/Dg-branes describe a certain 
kind of phase transitions in gauge theories [3], and its motion corresponds to the motion of quarks, 
mesons and baryons in the quark-gluon plasma [3] . It is also an interesting observation that the brane- 
black-hole system captures several aspects of the 'merger' of black holes appearing in the black-hole 
black-string system in Kaluza-Klein spaces (5l[6]. Besides these applications, the minimal surfaces, i.e., 
the minimizer of the Dirac-Nambu-Goto action, is used as a technical tool to prove some important 
theorems in general relativity such as the positive mass theorem [7]. 

In every context of the above applications of brane-black-hole system, a natural question arises as 
to which equilibrium state is most stable or what is the final state of dynamics such as gravitational 
scattering and capturing processes. The answer, of course, will depend on both the symmetries of the 
background black-hole spacetime and the boundary conditions of the brane at the asymptotic regions 
(infinity and horizons). As a specific example, let us consider a spherically symmetric static black hole 
as a background geometry. A naive expectation is that the final state of the dynamics of an isolated 
brane-black-hole system is the configuration in which the brane is trapped on the equatorial plane 
of the black hole, because the equatorial plane, which has the highest symmetries, will always be an 
extremizer of the action and attractive force of the black hole seems to put such a brane at the bottom 
of a concave potential. 

As we will see soon, however, the problem turns out to be more subtle practically than one expects 
as above. To assess the linear stability of the brane at the equatorial plane, one is forced to evaluate 
a one-dimensional Schrodinger equation as occurred in a perturbative stability analysis of black hole. 
As known well, a careful treatment is required to conclude the stability even for higher dimensional 
static black holes possessing maximally symmetric horizons [3 [9]. In particular, the presence of a 
cosmological constant and a charge/momentum makes the potentials in the Schrodinger equation 
non-trivial and complexity increases as dimensionality becomes higher. The intricacy of the effective 
potentials is the major obstacle for the analytic study of stability. 

In addition, we should remind some facts known in the mathematical studies of minimal surfaces 
in flat Euclidean space E"", which has a long history since the pioneering works by Young and Laplace. 
The plane in E'^, which may be regarded as a counterpart of the brane on the equatorial plane in 
our setting, is the only globally well defined single valued graph of the minimal surface equation, 
which is known as the Bernstein theorem, and its expected extension to general dimension is known as 
the Bernstein conjecture. Although the spacetime dimensionality appears to play no role to prevent 
this conjecture, this is not the case [lOj . The minimal surface need not be smooth necessarily in 
arbitrary dimensions, and there exists a rather explicit example called minimal cone for n > 8, which 
is curved and even singular at the apex. Quite recently, it was pointed out that this failure of the 
Bernstein conjecture and the existence of minimal cones had involved a variety of gravitational and 
non-gravitational systems [11], including the stable (local) cone solution of brane in the brane-black- 
hole system [6]. 

As far as the authors know, an initial attempt to address the stability of brane trapped on the 
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equatorial plane of a black hole was made by Higaki et al. [12], and this appears to be the only 
study in the literature. They discussed the 4-dimensional Reissner-Nordstrom-de Sitter black and 
found a conclusive result that the electromagnetic charge and positive cosmological constant tend 
to destabilize the brane. As mentioned above, however, the stability of minimal surfaces/branes in 
general dimensions has non-trivial features, and the results of perturbation analysis would shed light 
on this issue. Motivated by this and recent growing interest in multidimensional dynamics of gravity, 
this paper explores the stability of branes trapped on their equatorial planes of d-dimensional black 
holes {d > 4), including the case of asymptotically anti-de Sitter black hole while switching off the 
electromagnetic field for simplicity. 

The organization of this paper is as follows. In the next section, we briefly review the covariant 
perturbation method of a test brane. In Sec. [3] we specify the background spacetime and reduce the 
covariant equation to a single eigenvalue problem in radial direction. In Sec. H] the (in-)stability of 
brane is examined. The final section is devoted to discussions. In Appendix El a prescription to give 
the boundary condition of perturbations is given. We use the geometrical unit c = Gd = 1, where Gd 
is the d-dimensional gravitational constant, and (—,+,+,...) sign convention. Our notation of the 
cosmological constant is that the Einstein equation in (i-dimension without a matter content takes the 
form of Rfj,u = 2Ag^,^/(d - 2), and we rewrite A = €{d - 2){d - l)/(25^) where e = 0, ±1 and h (> 0) 
is a length scale. 



2 Covariant Perturbation Method 

Let us consider a timelike hypersurface embedded in the background spacetime (M, g^y) as x'^ = 
X^^{^°-), where {fj,,u = 0, 1, . . . , d — 1) and (a, 6 = 0, 1, . . . , d — 2) are the coordinates of the 
spacetime and the hypersurface, respectively. The Dirac-Nambu-Goto action is 

I=-r j d'^-^i^, (1) 

where r is the tension of the brane and 7 is the determinant of induced metric given by 

lab := gf.udaX^'dbX'' . (2) 

The equation of motion derived from the action is expressed as the vanishing trace of extrinsic curva- 
ture, 

r'Kab = 0. (3) 

Specifically, the brane motion is described by minimal surfaces. 

The perturbation of a background brane equilibrium can be thought of as a test scalar field on 
the background brane regarded as a (d— 1) -dimensional submanifold (Y^,-fab)- A systematic covariant 
treatment was given in [13]. The quantitative measure of brane bending is expressed by projecting 
the variation of the brane 5X^ into normal direction of the background equilibrium, 

^> := 71^,6 X^' . (4) 
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Varying the equation of motion ([3]), it turns out that the hnear perturbation equations simphfy to a 
scalar field equation on the hyper surf ace, 

n(T)$ - - h^'^R^^'^ $ = , (5) 

where /i^jy := g^v—n^riy is the projection tensor onto the worldsheet. n^"^) = (— 7)~^''^9a[(— t)^''^ T^^c^fc] 
and are the d'Alembertian operator and the Ricci tensor associated with the induced metric 7^^, 
respectively. 



3 Reduction to an Eigenvalue Problem 

As a background spacetime, we consider the static, spherically symmetric d-dimensional spacetimes 
{d > 4) whose line element is 

gf.^dx'^dx'' = -/(r)df2 + f{r)-^dr^ + {dO^ + sin^ ^dfi^.g) , (6) 

where dO^_3 is the line element of a unit (d— 3)-sphere. Let us consider the Dirac-Nambu-Goto {d— 1)- 
dimensional brane {T,,^ab) with an 0{d — 2) invariance embedded in the background ([6]). Denoting 
the brane configuration by 9 = 9{r), the induced metric reads 

labdede = -fir)dt^ + [fir)-' + r^e'\r)] dr^ + W 9 d^j^, . (7) 

Without writing down the explicit form of Dirac-Nambu-Goto equation of motion, one might deduce 
that the equatorial plane (9 = it/2) always extremalizes the action [since the equatorial plane is 
geodesic in the (r, 0)-plane]. This can be easily checked as follows. For the present brane configuration, 
the Dirac-Nambu-Goto action computes to give 

/ = -rOrf_3(t2 - h) j dr {rsm9)'''^V^ + fr^0'\ (8) 

where ^d~3 is the volume of a unit {d — 3)-sphere, and {t2 — ti) is a time interval. Then, varying this 
action, one obtains an equation of motion, 

r + ^ [2(. - 2)/ + ./'] S" - i^,-' + + -1^ . (.) 

Thus, 9{r) = Tr/2 obviously solves this equation, as we desired to show. Hereafter, we shall concentrate 
on the brane located at the equatorial plane and consider its linear perturbations. 

The variables in equation ([5]) are separable by setting <1> = r~('^~^)/^x(r)e~*'^*V(r2) and introducing 
the tortoise coordinate dr^, = dr//. Here, Y{Q) is a spherical harmonic on a (d — 3)-sphere (see, 
e.g., [8]) satisfying 

A(«f-3) + + d - 4)1 y = , £ = 0,1,2,..., (10) 
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where A*^"^ is the Laplace-Beltrami operator on the unit {d — 3)-sphere and i is the multipole index. 
It then follows that the wave equation reduces to the Schrodinger form, 



dr^ 



+ V{r)x = u;'x, (11) 



with an effective potential, 
^(0 = ^ 



£(£ + d - 4) - (d - 3) + i^^lf^/ + 



(12) 



We have hitherto nowhere used Einstein's equations to derive the perturbation equations. From 
now on, to make the discussion reasonably focused, we specialize to the cases where the background 
spacetime satisfies the vacuum Einstein equation with or without a cosmological constant, R^i, = 
e{d— l)g^,^/b'^, where e = 0, ±1 and b (> 0) is a curvature radius. That is, we consider the case where 
/(r) takes the form, 

where rg (> 0) is a constant corresponding to a mass of the black hole. The case of e = and e = 
+1(— 1), respectively, corresponds to the d-dimensional Schwarzschild black hole and Schwarzschild- 
dS(AdS) black hole|l| For the simplest Schwarzschild case, we first remark that the range of r^: is 
complete — oo < < oo and hence the outside wedge of a black hole is globally hyperbolic, indicating 
that the dynamics is uniquely determined by the initial data. The parallel will be inferred for the 
Schwarzschild-dS case if we recognize r^, — > oo at the cosmological horizon rc instead of infinity, where 
Vc is the larger root of /(r) = 0. In contrast to these two cases, the range of r+ is incomplete, which 
we take — oo < < 0, for the Schwarzschild- AdS, implying that the outside region of a black hole 
ceases to be globally hyperbolic. This distinguished property demands a special care for the dynamics 
of X in the e = — 1 case. 

For the metric (jlSp . the potential is written explicitly as 



M{i + d - 4) + {d - 5){d - 3) + {d - 3){d - 1) (^) - e{d' - 1) ( ^ 



(14) 



Inspecting (jl4p . one finds a universal property that the potential V is bounded below and vanishes 
at the event horizon —00, independent of the sign of the cosmological constant. However, 

behaviors of the potential V away from the event horizon is very sensitive to e. For e > 0, the 
potential vanishes at r* ^ 00 (infinity for e = and cosmological horizon for e = +1). A distinct 
feature of a positive cosmological constant is that it will give a negative contribution to the potential 
[the last term in equation ()14p ]. By virtue of this, V approaches to zero from below as — > 00 and 
causes an instability, as we will show later. Whereas, the potential V diverges at infinity in the case 
of e = — 1. This illustrates that the negative cosmological constant acts as a 'confining box.' In either 
case, the potential V fails to be positive-definite in general. 

In the rest of this paper, we examine the stability of the brane for each background. 



^ More precisely, the metric function (|13p for e = 1 describes a Schwarzschild-dS black hole provided ro satisfies 
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4 Stability Analysis 



In order to conclude the linear stability, we must show that the Schrodinger equation Ax = uP'x 
[equation (jlip ] admits no normalizable negative mode solutions uP' < 0, where we have defined A := 
~(P /dr'^ + V . The operator A is elliptic and identified as a Hamiltonian operator on a Hilbert space of 
square-integrable functions L^(r*,dr*) on a static timeslice. To establish the stability, it is necessary 
to show that A is a positive, self-adjoint operator with L^(r*, dr*)-norm. 

Before embarking on the stability analysis, let us digress here and first focus on the dynamics of a 
scalar field x in this background. Since the outside wedge is globally hyperbolic for the Schwarzschild(- 
dS) spacetime, the dynamics of x is well-posed in this region. Since the potential V vanishes as 

— > iboo, the only possible boundary condition to obtain a normalizable solution is the Dirichlet 
boundary condition x = at — > ±00 (see Appendix |A]) . Accordingly, the self-adjoint extension Ae 
of A is unique, that is to say, the Hamiltonian A is essentially self-adjoint. Hence all that remains to 
solve for the stability is to show the positivity of Ae- This is a standard prescription to pursue the 
dynamics in globally hyperbolic spacetimes. 

For the Schwarzschild-AdS case, on the other hand, the domain of outer communication is no 
longer globally hyperbolic. This means that the ordinary Cauchy evolution determines a solution of 
evolution equations only within the domain of influence for a given initial data slice. Nevertheless, this 
difficulty can be cured for a static spacetime and it is always possible to find a sensible dynamics-and 
this is essentially the unique recipe for defining dynamics under quite reasonable conditions-beyond 
the domain of dependence of initial data slice [15]. To define a unitary dynamics throughout the 
non-globally hyperbolic static spacetime, we need a self-adjoint extension of A, which may or may not 
be unique. As is well known, choosing a self-adjoint extension is equivalent to choosing a boundary 
condition. 

Let us determine possible boundary conditions for e = — 1. Since the potential vanishes at the 
horizon, the normalizability singles out a unique boundary condition therein. Henceforth, our primary 
concern is the boundary condition at infinity. Near infinity r^: ~ —b'^/r ~ 0, ^ behaves as 

this yields the asymptotic solution, 

X-cM'^'^^' + ---) + C.{r:'''-'^/' + ...). (16) 

The normalizability requires C2 = 0, viz, V is in the limit point case at infinity. Therefore, the 
self-adjoint extension of A is found to be unique even for e = — 1. This observation is central to the 
stability discussion developed below. 

Given a self-adjoint extension Ae of A, let us next move on to the issue on the positivity of Ae- 
If the potential V happens to be positive, Ae corresponds to a Friedrichs extension and the stability 
against linear perturbation immediately follows [H]. As we have seen, however, the potential function 
V, equation p^ . is not positive-definite, so it is far from obvious whether Ae has positive eigenvalues. 
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To see the positivity of Ae, it is of great benefit to employ the idea proposed by Ishibashi and 
Kodama [9]. Fohowing their algorithm, the potential in the Schrodinger equation can be 'deformed' 
in such a way that a newly deformed potential is positive-definite. Let S denote an arbitrary function 
of r and 



D 



_d_ 

dr* 



S 



be a derivative operator. Straightforward calculation shows that 



(X, Ax) 



L2 



XDX 



boundary 



+ / dr* 



Dx 



and 



Vs = V-f—S-S^ 
ar 



(17) 



(18) 



(19) 



It deserves to remark that this expression is formal because there is in general no guarantee that the 
boundary term and the integrand are both finite. If one can show that the boundary term vanishes 
and if one can find an appropriate S that makes Vs positive, the Hamiltonian operator A turns out 
to be a positive operator, which has at least one positive self-adjoint extension corresponding to the 
Friedrichs extension. Utilizing this formula, we are thereby able to demonstrate the positivity oi Ae, 
i.e., the stability of test brane perturbations. 



4.1 Schwarzschild(-AdS) Background 

Let us begin by the analysis of e < 0. As we have shown that the boundary condition is necessarily of 
Dirichlet-type, it then follows that the boundary term in equation (jl8p indeed drops off and > 
implies a stable dynamics. We can accomplish this by the choice 



d-3 
2r 



giving rise to a positive potential for any ^ > 0, 



(20) 



(21) 



Then, the test brane perturbation in Schwarzschild(-AdS) black hole is stable for any d > 4 dimen- 



sionsJ 



4.2 Schwarzschild-dS Background 



According to [12], the four-dimensional Schwarzschild-dS black hole exhibits an instability only for 
the i = mode. We now show that this result continues to be true for all d > 4. 



^One can show that the string motion in the background of BTZ black hole {d — 3, e = —1) is also stable. 
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As the domain surrounded by event and cosmological horizons is globally hyperbolic, equation (fTHl 
receives no boundary contribution. We find that 



leads to the deformed potential 



VsJi^M±izM. (23) 



So the system is stable in d > 4 dimensions against for all ^ > 1 modes. A possible unstable mode 
exists only for the s-wave|f] 

We will next show numerically that the s-wave actually excites an unstable mode. To perform 
numerics, we have to make dimensionless quantities. Here, the normalization by ro seems not so 
useful because it is a measure of the mass, which has the dimensional dependence as Tq"^ ~ G^Af . To 
compare black holes in different spacetime dimensions, it is more advantageous, instead of ro and 6, 
to use Tc and r/j, which are, respectively, the radii of the cosmological horizon and the event horizon. 
Simple calculations show that they are related by 



d-i ' b - ,^ Jj_-^ . (24) 



One may recognize immediately that the background becomes the Schwarzschild black hole as the 
ratio rc/vfi tends to infinity. 

We show the values of growth rate a := Im(c<j) of perturbations in typical dimensions in Table [ifl 
When the ratio rc/r^ is fixed, the growth rate ar^ becomes higher as the dimension gets higher. The 
brane on the equatorial plane of the higher dimensional Schwarzschild-de Sitter black hole is more 
unstable than that of the lower dimensional one with the identical value of ratio Vc/rh. It can be 
understood by the form of potential in each dimension. When the spacetime dimension is fixed, the 
growth rate arh becomes lower as the ratio Vc/rh gets bigger. In other words, smaller cosmological 
constant reduces the instability of the system. This tendency is compatible with the fact that the 
system in the case of the Schwarzschild black hole background is stable. 

To sum up, we found that the positive cosmological constant strengthens the instability of branes 
on the equatorial plane for the arbitrary dimensional Schwarzschild-de Sitter black holes. We can also 
say that the higher dimensional background has a destabilization effect of the system. 



^It is worthwhile to emphasize that equation (|23|) can be derived from equation (|12|) without any reference to the 
particular form of /. Thus, all the spherically symmetric static black-hole type backgrounds of the form ([6| admit a 
stable brane dynamics for 1>1 modes. 

*The proof of uj being pure imaginary is relegated to Appendix. The numerical results presented in Table [T] were 
obtained by two independent shooting calculations. However, the numerical credibility becomes worse as the dimension- 
ality increases, d > 10, in particular for large rc/rn. It would be interesting to apply the matched asymptotic expansion 
method for the large-d limit [17], by which one may be able to determine the growth rate analytically. 
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Table 1: Numerical values of the growth rate a = Im(a;) of the i = perturbation of brane in the 
Schwarzschild-de Sitter background. The horizon radius is set to unity. 





d = 4 


d = 5 


d = 6 


d = 10 


d = 11 


2 


1.332 X 10"^ 


2.429 X 10"^ 


3.277 X 10^^ 


4.762 X 10"^ 


4.865 X 10-1 


20 


4.435 X 10"^ 


4.957 X 10"^ 


4.997 X 10"^ 


5.00 X 10"^ 


5.00 X 10"^ 


100 


9.750 X 10"^ 


9.996 X 10"=^ 


1.000 X IQ-^ 


1.00 X 10-^ 


1.00 X 10^^ 



5 Discussions 

We have seen that the branes rest at the equatorial plane of the (i-dimensional {d > 4) Schwarzschild 
and Schwarzschild-AdS black holes are stable for the linear perturbations, while for the Schwarzschild- 
dS background the brane is unstable against the s-wave perturbation. The covariant perturbation 
equation of brane is translated into the the eigenvalue problem in the Schrodinger form after the 
separation of variables. Then, the stabilities were proved by showing the positivity of the Hamiltonian 
operator by means of a deformation technique. In the asymptotically de Sitter case, the positivity of 
Hamiltonian operator only for the non-spherical perturbations {I > 1) was shown, which means that 
one has to solve the eigenvalue problem explicitly to pursue the behavior of spherical mode. To this 
end, we numerically solved the equation for the spherical perturbation and showed that the mode is an 
unstable one in every dimension, growing exponentially in time. One can recognize this coming from 
the fact that large £ has the stabilizing effect. The dangerous modes originate from the lower angular 
eigenvalues since the effective potential becomes higher as £ increases. This result is reminiscent of 
the fact that the Gregory-Laflamme instability of black branes exists only in the s-wave sector |18] . 

For definiteness of our argument, we have limited our consideration to a neutral black hole as a 
background black hole. The inclusion of electromagnetic charge is straightforward. The electromag- 
netic charge will leave a negative imprint on the effective potential (|12|) as in an analogous fashion 
in four dimensions [12], where an instability was explicitly shown. Thus there seems no compelling 
reason to stabilize the brane motion in the charged case. As we said, this instability occurs only for 
the s-wave in any spacetimes of the form ([6]) such as a black hole in different gravitational theories. 

In this paper, we have ignored the backreaction of brane self-gravity and the thickness of brane, 
both of which cannot be negligible in general situations. If one takes into account the backreac- 
tion, it is known that the domain wall exhibits repulsive gravitational fields [19]. Regarding the 
thick branes interplaying with black holes, several equilibrium solutions were constructed numerically 
or analytically [20\ I21j . The introduction of the thickness is known to result in the expulsion of 
branes/domain walls by an extremal/near-extremal black hole. Thus, it would be interesting to see 
how the competition between the attractive and repulsive forces appearing in more realistic modelings 
of brane-black-hole system affects the stability of the brane examined in this paper. 

While we have restricted ourselves to the linear stability of brane and found, for example, the 
stability in the Schwarzschild background, it is known that a certain sort of non-linear perturbations 
to the Dirac-Nambu-Goto brane in the Schwarzschild background, which mimic the recoil of black hole 
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due to the Hawking emission, evolves and results in a reconnection of the bent brane, corresponding to 
the eventual escape of black hole into the bulk [22] . It is noted that our results are not inconsistent with 
those in [22] since the asymptotic boundary conditions are different. In a similar vein, the existence of 
a critical escape velocity of the black hole from the brane was investigated in [23]. Interestingly, the 
existence of critical escape velocity depends on the number of codimensions. Although we have been 
concerned only with the codimension-one branes in this paper, it would be interesting to investigate 
how the number of codimensions affects the stability of the brane, which can be examined with the 
covariant perturbation methods developed in |24j . 
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A Boundary Conditions 



Here, we assemble several basic facts about the behaviors of physical perturbations in the case of 
e > 0, which are used in this paper and might be helpful for readers. 

For the Schwarzschild case (e = 0), — > — oo as r — > r/j (r = r/j is the event horizon) and 
— > +00 as r — > +oo. For the Schwarzschild-dS case (e = +1), — > — oo as r — > r/j again and 
r* — > +00 as r — > Tc (r = Tc is the cosmological horizon). One can see that the potential (jl4p vanishes 
as r* — > iboo. Thus, the linearly independent asymptotic solutions are the plane waves, e^*"^^'*. Since 
we assume the time dependence of x oc e"*"^*, two modes e~*'^^* and e"^''* correspond to the ingoing 
and outgoing waves, respectively. Hence, we impose the following boundary conditions on a physical 
mode, X ~^ e^*'^''* as r^: — > ±oo. 

As far as an unstable mode, corresponding to Im(Li;) > 0, is concerned, we can show that u> is pure 
imaginary as follows. Multiplying the complex conjugate of X; with equation ([TT]) . and integrating 
it, we have 



(x,^x)l2 



dx 



dr* 



+ V\x\^ dr. 



X 



dx 
dr* 



u> 



Ixl dr*. 



(25) 



Since the surface term in equation (j25p vanishes due to the exponential suppression at infinity, the all 
integrals in equation (|25p are real. Thus, must be real, and furthermore uj must be pure imaginary 
[since we are looking for the unstable mode, lm{uj) > 0]. Therefore, we can write as uj = ia (a > 0) 
and equation ([25|l reads 

2 



+ 00 



dx 



dr* 



+ V\x 




1x1 dn 



(26) 



What we can see from this equation is that if the potential V{r) is positive definite, the right hand 
side is positive and equation (|26p makes no sense, implying the absence of unstable mode. In terms 
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of (T, the asymptotic boundary conditions mentioned above are rewritten as x e as r^, — > =Fc>o, 
which are used in the numerical analysis. 
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